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ABSTRACT 

FS CMa type stars are a recently described group of objects with the B[e] phenomenon that exhibit strong 
emission-Hne spectra and strong IR excesses. In this paper we report the first attempt for a detailed model- 
ing of IRAS 00470-1-6429, for which we have the best set of observations. Our modeling is based on two key 
assumptions: the star has a main-sequence luminosity for its spectral type (B2) and the circumstellar enve- 
lope is bimodal, composed of a slowly outflowing disk-like wind and a fast polar wind. Both outflows are 
assumed to be purely radial. We adopt a novel approach to describe the dust formation site in the wind that 
employs timescale arguments for grain condensation and a self-consistent solution for the dust destruction sur- 
face. With the above assumptions we were able to reproduce satisfactorily many observational properties of 
IRAS 00470-1-6429, including the H I line profiles and the overall shape of the spectral energy distribution. Our 
adopted recipe for dust formation proved successful in reproducing the correct amount of dust formed in the 
circumstellar envelope. Possible shortcomings of our model, as well as suggestions for future improvements, 
are discussed. 

Subject headings: stars: emission-line — stars: early-type — stars: - circumstellar matter — stars: individual 
(IRAS 00470H-6429) 
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1. INTRODUCTION 

This series of papers is devoted to studying objects from the 
Galactic FS CMa type group. The group comprises about 40 
objects with the B[e] phenomenon, which refers to the simul- 
taneous presence of forb idden lines in the spectra and strong 
exces ses of IR radiation (lMiroshnichenkol2007t iLamers et al.l 
119981) . Its members show properties of B-type stars with typ- 
ical main-sequence luminosities [log(L/L0) ~ 2.5-4.5]. On 
the other hand, they exhibit extremely strong emission-line 
spectra that are on average more than an order of magnitude 
stronger than those of Be stars of similar spectral types. The 
group includes objects that were discovered long ago, but 
their nature and evolutionary state wa s not satisfactorily de- 
termi ned, as well as newly found ones (iMiroshnichenko et alj 
l2007h . 

One of the main problems in understanding those objects 
has to do with the presence of significant amounts of hot dust 
in their circumstellar region that cannot be left from the time 
of the stars' formation. It is also hard to explain the dust 
formation using typical mass loss rates from single stars of 
an appropriate mass range (~ 5-2OM0). Our observations 
show that over 30 per cent of the group objects show signs 
of secondary companions in the spectra (e.g., Li I 6708 A 
absorption line). Therefore, one might assume that the large 
amounts of circumstellar matter near FS CMa type objects 
may be a consequence of the binary evolution, although direct 
evidence of the mass transfer between the stellar companions 
is not observed. 

|carciofi@usp.br] 



The main goals of this series of papers are to obtain high 
quality observational data for the group objects, to describe 
the objects' properties in detail, and to model the latter in or- 
der to put tight constraints on parameters of the stellar com- 
panions and those of the circumstellar material. The ob- 
served properties of the first object for which we collected 
a good wealth of data, IRAS 0470-H6429, were presented 
by IMiroshnichenko et al.l (120091 hereafter Paper III). In the 
present paper we make use of the Monte Carlo code HDUST 
to determine the physical parameters of this intricate system. 

2. MODEL DESCRIPTION 

From the discussion of the available data in Paper III, it be- 
comes clear that IRAS 00470-1-6429 is a very complex object, 
with a large and dense CS envelope composed of both ionized 
gas and dust. It is expected, therefore, that any model that 
successfully reproduces the observed properties must also be 
complex. 

This first attempt of modeUng of IRAS 00470H-6429 
was conducted usi ng th e computer code HDUST 
JCarciofi & BiorkmanI EOOSl l2006t ICarciofi et al.] l2008l 
120061 1200411 HDUST is a three-dimensional Monte Carlo 
radiative transfer code that combines the full non-local ther- 
modynamic equilibrium (NLTE) treatment of the radiative 
transfer in gaseous media with a very general treatment for 
CS dust grains. This code has been already successfully 
applied to study the CS d isks of Be stars (see, for instance, 
ICarciofi et al]l2007ll2009l) . 

In HDUST, the CS environment can be represented by two 
components: a gaseous region, containing hydrogen and free 
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TABLE 1 

Model Parameters for IRAS 00470+6429 



Parameter 


Value 


Type 


Stellar- Parameters 


Spectral Type 


B2V 


fixed 


R 


6R0 


fixed 


ne 


20000 K 


fixed 


L 


5 100 Lq 


fixed 


Geometrical Parameters 


(■ 


84° 


free 


d 


llOOitlOOpc 


free 


Wind Parameters 


dM(90°)/cm 


(1.4— 1.6) X lO^'^Mgyr"' 


' sr ' calculated 


dM{Q)/dn 


< 3 X 10-' Mq yr-l sr*' 


free 


A, 


> 19 


free 


Voo(90°) 


24 km S-' 


calculated 


Voo(O) 


1200 km S-' 


fixed 


A2 


-0.98 


free 


/9(0) 


4 


free 


A3 





free 


m 


92 


free 


Rw 


> 10-'^ «0 


fixed 


Common Dust Parameters 


Composition 


Amorphous silicate^ 


fixed 


n 


-3.5 


fixed 


Gas-to-dust ratio 


200 


fixed 


Dust Model 1 


^'min ^Imux 


0.05 — 0.25 pm 


fixed 


Pdust 


1.5gcm-3 


free 


-'destruction 


1 500 K 


free 


'^destruction 


540 R 


calculated 


V 


0.02 


free 


E{B-V)is 


1.1 


free 


Dust Model 2 


^min ^max 


0.05 — 0.10 ^m 


fixed 


Pdust 


0.15gcm-3 


free 


-'destruction 


1 500 K 


free 


'^destruction 


540 R 


calculated 


V 


0.02 


free 


E(B -V)is 


1.1 


free 


Dust Model 3 


^min ^mux 


1— 50 pm 


fixed 


Pdust 


0.1 gcm-3 


free 


-'destruction 


1600K 


free 


'^destruction 


79 R 


calculated 


V 


0.02 


free 


E(B-V)is 


1.2 


free 


'The optical constants 


i oflOssenkoDt. Hennine & MathisI 09921) were us 



electrons, and a dusty region. The solution to the problem 
is obtained in two steps. The first step consists of an itera- 
tive scheme to solve the coupled problems of radiative trans- 
fer, radiative equilibrium and statistical equilibrium. For the 
gaseous shell, HDUST calculates the NLTE hydrogen level 
populations and the electron temperature as a function of po- 
sition, taking into full account all the relevant physics (con- 
tinuum and line transfer, collisional processes, etc.). For the 
dusty shells, HDUST solves the radiative equilibrium prob- 
lem, thereby obtaining the equilibrium temperature of the dust 
grains a s a function of position. The interested reader is re- 
ferred to lCarciofi & BjorkmanI (120061) for details of the Monte 
Carlo NLTE solution and ICarciofi et al.l ( 120041) for informa- 
tion about the dust implementation. Once the first step is com- 
pleted and the state variables of the CS material are known, 
HDUST is run in a post-processing mode to calculate the emer- 



gent flux (spectral energy distribution, line profiles, synthetic 
images, etc.) for arbitrary lines of sight. 

A few words about an important feature of the code are war- 
ranted. For a given gas pl us d ust distribution, as, for instance, 
the one adopted in Sect. 12.11 HDUST solves for the gas and 
dust properties simultaneously and self-consistently. Thus, 
important effects, such as the shielding of the dust grains from 
the stellar ionizing UV photons by the CS gas, are included 
in the calculations. This feature allowed us to determine for 
each model the dust destruction surface, which is the latitude- 
dependent distance from star beyond which the dust equilib- 
rium temperature is below an adopted destruction tempera- 
ture. This feature of the code proved to be a key factor for the 
model we describe below. 



2.1. A Parametric Model for the Circumstellar Envelope 

For the problem at hand, HDUST must be provided with a 
description of the properties of the central star (effective tem- 
perature, Teff, and radius, R) and the physical properties of the 
CS envelope (density and velocity distribution of the gas and 
the characteristics of the CS dust, such as composition and 
grain size distribution). In addition, two geometrical param- 
eters must be specified: the distance to the star, d, and the 
inclination angle, /, which is the angle between the observer's 
line of sight and the symmetry axis of the system. Finally, 
the interstellar reddening, E{B-V)i^, must also be specified. 
Even though it is not an intrinsic parameter of the system, we 
include it as a model parameter because of the uncertainty in 
its determination, mainly owing to the difficulty to separate 
the interstellar from the circumstellar reddening. 

A long standing issue regarding the stars with the B[e] phe- 
nomenon is that very little is known about the proc esses that 
lead to the formation of their CS environment (see IZickgran 
I2OO6I for a review). Also, in spite of the effort of numer- 
ous observers and theoreticians, the most fundamental prop- 
erties of the CS envelope remain quite uncertain. In the 
case of supergiants B[e] (s gB[e1), a subclass of the stars 
with the B[e] phenomenon (iLamers et al.lll998l) . it is gener- 
ally believed that their CS envelope is bimodal, composed 
of a slow ly outflowing, possib ly rotating disk and a fast po- 
lar wind dZickgraf et al.lll985h . Two scenarios invoked to ex- 
plain this bimodality are th e rotationally induced bi-stability 
(ILamers & P auldrach"1991l) and rotationally supported Kep- 
lerian outflows (Porter 2003). 

CS envelopes of FS CMa objects are even less understood, 
since very little theoretical work has been done with those 
objects so far For this initial theoretical effort we adopt a 
model for IRAS 00470-1-6429 that has two key assumptions. 
First, we assume that the star has a luminosity similar to that 
of a main-sequence star with the same spectral type, i.e., the 
stellar luminosity is in the range of a few thousand solar lu- 
minosities. We believe this assumption is well supported by 
the available data (Paper III). Second, we assume that some 
sort of bimodal, two-component wind is present, similar in 
structure to the CS envelopes of sgB[e]. The results shown in 
Sect. ?? largely corroborate the above assumptions, but pos- 
sible shortcomings are discussed in Sect. |4] 

We adopt the following parametric form for the CS density 
and velocity structure. We assume that the underlying phys- 
ical reason for the bimodal envelope is that the mass loss is, 
somehow, enhanced around the equator We express the stel- 
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lar mass loss rate per unit solid angle as 



dMie) dMiO) 



[l+Aisin^C^)] , 



(1) 



dn dn 

where 6 is the colatitude measured from the pole. Ai is a 
parameter that controls the ratio between the equatorial and 
polar mass loss rates 



dM{90°) dn 

Al = : ■ 1 

dM(0)/dn 



(2) 



Since the mass loss rate is enhanced at the equator, Ai > 0. 
The parameter m controls how fast the mass loss rate drops 
from the equator to the pole. Defining the disk opening angle, 
A^op, as the latitude for which dM(d)/dft has dropped to half 
of its equatorial value, we have 



A^op = sin 



■1 



2A, 



l/m' 



(3) 



For the velocity structure, we assume that the gas motion is 
purely radial, with a latitude-dependent radial velocity given 
by 

Vr{r,0) = VQ + [voo(e)-vo](l-R/rf^'\ (4) 

where r is the distance to the center of the star. For a given 
^, eg. (Bli is a standard /3- law for radiatively driven winds 
jLamers & Cassinenilll999t) . We assume that both the wind 
terminal velocity, Voo, and the acceleration parameter, /3, are 
functions of 9, and have the same latitudinal dependence as 
the mass loss. Thus, 



and 



VooW = Voo(O)[l+A2sin'«(0)] , 
/3(6») = /3(0)[l+A3sin'"(6i)] . 



(5) 



(6) 



As before, A2 and A3 control the ratio between the values of 
each quantity at the equator and at the pole. 



and 



M = 1 , 

Voo(O) 

/3(90°) 



A,= 



/3(0) 



1 . 



(7) 



(8) 



Since the terminal velocity at the equator is smaller than at the 
pole, A2 < 0. 
Assuming radial outflow, the mass continuity equation is. 



dMiO) 2 



dn 



r'p(r,e)v{r,e), 



so we can write the CS gas mass density distribution as 

dM(9)/dn 



P(r,e): 



rMr,0) 



(9) 



(10) 



2.2. Dust Condensation Site 

Having chosen a prescription for the CS gas, we must now 
provide a site for dust formation and describe the properties 
of the dust grains that are for med. In an in itial study of the 
continuum emission of sgB[e1. lPorteil ( 120031) used two criteria 
to determine the dust formation site. In his models, dust was 
formed where the gas temperature, which he assumed to be an 
ad hoc function of the distance from the star, was lower than 



a given value (w 1 500 K) and the gas density was larger than 
a critical value. 

In our models, the coolest part of the envelope is in the 
denser equatorial regions, an expecte d result based on pre- 
vious works (e^. 'Zsar go et al.l 120081 : ICarciofi & BjorkmanI 
I2OO6I see Fig. |5l below). However, even in the midplane 
the gas temperatures never fall below about 4000 K. This 
temperature is much larger than the temperature needed for 
dust grain condensa tion to occur, which is thought to be be- 
low about 1500 K (lPortedl2003h . It is likely that the gas 
temperatures predicted by HDUST are overestimated, because 
no cooling by heliu m or other metals is included. However, 
IZsargo et al.l (l2008h . in a study of the properties of the CS en- 
velope of sgB [e] , found that the temperature never falls below 
2000 — 3000 K, even though they used more realistic cool- 
ing terms. These authors argue that a possible mechanism for 
further cooling the gas is adiabatic cooling, or an even more 
complex chemical composition. 

Because of the above, it is clear that the gas temperature 
cannot be used in our models to trace the dust formation site 
and we use, for this p urpose, a different approach than the one 
used by lPorteii(l2003h . We define the dust formation site using 
two complementary criteria. If dust is to be present at a given 
location of the CS envelope, 

i) The equilibrium temperature of the dust grains (not the 
gas temperature) must be smaller than a given grain destruc- 
tion temperature, rdestmction- 

ii) Furthermore, the gas density must be larger than a criti- 
cal value, such that grain growth can occur 

The first criterion was implemented in HDUST using a very 
simple procedure. For a given model, we start from a guess of 
the distance from the star beyond which the dust grains begin 
to condensate. During the iterative process used for deter- 
mining the temperatures, whenever the dust temperature in a 
given cell becomes larger than Tdestmction, the code removes the 
dust from that particular cell. Conversely, if in a given dust- 
free cell the temperature of the dust, if present, would have 
been smaller than the destruction temperature, then that cell 
becomes a dusty cell. This simple procedure ensures a self- 
consistent determination of the location of the dust destruction 
radius as a function of latitude. Note that the dust destruction 
radius is highly model dependent: it is mainly controlled by 
the gas opacity in the inner part of the envelope and the dust 
model used. 

The second criterion for defining t he dust formation site fol- 
lows a simple argument outlined by iGail & Sedlmavil (119881) 
and adopted by Porter (2003). The idea is that a critical gas 
number density for the dust formation can be derived by com- 
paring two timescales, the timescale for chemical reactions, 
Teh, and the timescale for the gas expansio n, Texp. Grain for- 
mation and growth occurs if Teh < ^exp- iGail & Sedlmavd 
([198 8) found that this occurs if the number density of the 
growth species (silicon, in this case) is larger than the fol- 
lowing critical value 



nsiir, d) : 



v,-(r,0) 



10-"'n'rei(r,( 



(11) 



where Vrei is the relative velocity between the dust grain and 
an atom or molecule, which is of the order of the sound speed 
and, thus, position-dependent. Since «si of eq. ( fTTT ) was de- 
fined from simple physical arguments, and, therefore, is prob- 
ably just a rough estimate of the critical density, we introduce 
a scaling factor, 1], so that the condition for dust formation 
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becomes 






<f]nsi{r,9). 



(12) 



where me is the mass of the H atom and /i is the gas molecular 
weight, which is typically about 0.6 for an ionized gas with 
solar abundance, e is the relative abundance in number of Si 
to H, for which we adopt the solar value 3.6 x 10"^. The IR 
excess is a measure of the amount of dust that is formed in the 
CS envelope. Hence, the value of 77 can be found by matching 
the IR excess. The second criterion for the spatial distribution 
of the CS dust effectively confines the dust grains to the dense 
and slow equatorial regions. 

Finally, given a site for the dust in the CS envelope, we 
need to describe the properties of the dust grains, which re- 
quires several additional parameters. We assume that at each 
location of the envelope where dust can form, a certain frac- 
tion of the gas mass is converted into dust. This fraction, 
known as the gas-to-dust ratio, was assumed to be 200. We 
further assume that grains are spherical and their composi- 
tion is oxygen-based, i.e., the dust is comprised mainly of 
silicates. This assumption is supported by evolutionary argu- 
ments, since carbon is depleted due to the CNO processing cy- 
cle in massive stars, therefore ruling out carbon-based grains, 
such as graphite, for the object in question. Furthermore, IR 
spectra of nearly 20 FS CMa objects obtained with the Spitzer 
Space Observatory (IRAS00470H-6429 was not observed due 
to the presence o f a bright nearby IR source) show silicate 
emission features (iMiroshnichenko et al.ll2008l) . 

For this initial anal ysis we adopt a MR N law for the distri- 
bution of grain sizes (iMathis et al.lllQTTh . according to which 



f(a) = Ca", 



(13) 



where / is the fractional number of dust grains with ra- 
dius a and C is a normalization constant. In the standard 
MRN law, « = -3.5 and the distribution is confined between 
flmin = 0.05 fim and flmax = 0.25 (Um. In our models, we let 
both flmin and flmax to be free parameters in order to inves- 
tigate the possible existence of large dust grains (Paper III). 
Given a dust mass density and /(a), we further need the av- 
erage density of the grain material, pdust, in order to calculate 
the number density of dust grains. Since we do not know how 
fluffy the grains are, we let pdust be a free parameter. We note, 
however, that pdust and the gas-to-dust ratio are not indepen- 
dent parameters, since both are contained in the dust optical 
depth, the one quantity that effectively controls the radiative 
transfer in the CS dust shell. 

2.3. Adopted Grid Structure 

A few words are needed about the grid structure we adopted 
for the current problem. In HDUST, the CS envelope is divided 
into Nr radial cells, A^^ colatitude cells and N^ azimuth cells, 
and it is assumed that the state variables of the CS material 
(temperature and level populations) are constant in each cell. 

The density law of eq. (|4|i can present difficulties for a ra- 
diative transfer solver, because it has a large radial gradient 
close to the star and it may have a sharp enhancement at the 
equator, depending on the adopted value of ni. For the radial 
spacing we adopt cells sizes such that the radial electron scat- 
tering optical depth is the same for all cells, which results in 
very small cells close to the star and large cells in the outer 
envelope. This ensures that the inner regions of the envelope 
are sampled adequately. For the colatitude bins, we define 
the cell sizes is such a way that from the equator to the pole. 



along a ray with constant radius, the density at the center of 
consecutive cells falls in steps of2/N^. Therefore, when m is 
large, which corresponds to a disk with small opening angle, 
the cells are small at the equator and large at the pole. 

Two competing issues must be taken into account when 
choosing the values of A^^ and A^^. Since HDUST employs 
a statistical method for sa mpling cell-dependent quan tities 
such as radiative heating (see lCarciofi & Biorkmanl2006l) . the 
larger the number of cells the larger the computational effort 
necessary to reach a given sampling error in the rates. On the 
other hand, the grid structure must be fine enough to describe 
correctly the adopted density and to resolve possible regions 
of ionization changes, where a transition from ionized to neu- 
tral hydrogen occurs. The minimum number of cells neces- 
sary for a given problem can be estimated by running models 
with an increasing number of cells. If the cell spacing is prop- 
erly defined, there will be a certain number of cells beyond 
which no modification of the solution is found by further re- 
fining the grid. We found that A^^ = 100 (60 for the gas shell 
and 40 for the dusty shell) and A^^ = 40 are a good compro- 
mise between speed and accuracy. For the current problem 
the number of azimuth cells was set to 1 because axisymme- 
try was assumed. 

2.4. Modeling Procedure 

The adopted model described above has many parame- 
ters but for this initial study we chose to let just part of 
them be free parameters, while keeping the rest fixed. We 
chose as fixed parameters the ones which can either be suffi- 
ciently well-constrained by a model independent analysis, or 
are associated with the basic model assumptions defined in 
Sect. 12.11 Examples of fixed parameters are the stellar char- 
acteristics and the dust composition. 

Free parameters, on the other hand, are the ones about 
which no information can be readily obtained from the ob- 
servations and, therefore, require a detailed modeling. The 
main free parameters in our study describe the density and 
velocity of the CS gas and the latitude-dependent stellar mass 
loss rate. We believe that this approach is appropriate for this 
initial study, since it allows us to focus on the least known 
properties of the CS envelope. A list of model parameters for 
IRAS 00470H-6429 is given in Table [U 

2.5. Stellar Parameters 

The stellar parameters are all fixed in this study. As 
mentioned above, one of our key assumptions is that 
IRAS 00470H-6429 has a typical luminosity of a normal main- 
sequence B star. Following the analysis in Paper III we chose 
a spectral type of B2 V for the central star. This sets the radius 
to « 6 Rq, the effective temperature to k, 20000 K and the 
luminosity to « 5 100 L© (Harmanec 1988). We investigated 
the effects of changing %fi and L, and we found that varying 
Teff by ~ 10% (and thus L by w 40%) has only small effects 
on the results. For this initial study we considered a static, 
non-rotati ng central star. Thi s is not a limitation of HDUST 
(see, e.g. ICarciofi et al.l 120081) but, instead, reflects the fact 
that the observed line profile is controlled mostly by the CS 
gas (Sect. 13. it . 

For the stell ar spectrum we used the model atmosphere of 
iKuruczl (IT99I for Teft = 20000 K and log^ = 4E1 To describe 

' The use of a LTE model atmosphere is not a hmitation of the present 
work. More recent NLTE results could have been used, but the object is 
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the photospheric specific intensity as a f unction of dir ection 
we adopted the limb darkening models of lClareti ( 120001) . 

2.6. Geometrical Parameters 

The geometrical parameters (/, d and interstellar reddening) 
are all free parameters. Nothing is known about /, but our 
previous analysis allows us to impose some limits on d and 
E{B-V)is. Following Paper III, we adopt d in the range of 
1000 — 2000 pc and E{B -V)is in the range 0.8 — 1.2 mag. 
Given a model for the CS gas and dust, d and E{B-V)is are 
obtained by fitting the observed flux level and shape of the 
visible SED, respectively. 

3. RESULTS 
3.1. Circumstellar Gas 

Some information about the gas dynamics can be readily 
obtained by analyzing the emission line profiles. For instance, 
the broad blue absorption profile of the H I lines suggests that 
the absorbing material has large velocities in the line of sight 
and the blue edge of the absorption profile indicates the termi- 
nal velocity of the fast wind. We adopted for Voo(O) a value of 
1200 km s"' and set this as a fixed parameter in our modeling. 
We note that this value of the terminal velocity is consistent 
with what is expected for main-sequence B stars (iGrady et alJ 
[19871) . 

In addition, the relatively narrow peaks of the Balmer lines 
suggest that the bulk of the emitting material has small ve- 
locity along the line of sight. Supposing that most of the H I 
line emission comes from the equatorial region (an assump- 
tion corroborated by our detailed calculations) this suggests 
that the equatorial outflow velocities are very small compared 
to the velocities in the polar region. 

For the present work, hundreds of models were run cover- 
ing a wide range of model parameters. Before presenting our 
solution, it is useful to discuss some trends that became clear 
once a large range in the parameter space was covered. 



From our results we can safely rule out both large {Ado 



> 

15° orm < 15) and small (A6'op ^5° orm > 180) opening an- 
gles for the outflowing disk around IRAS 00470+6429. The 
opening angle, together with the equatorial mass loss rate, 
controls the strength of the H I line emission. For large open- 
ing angles, relatively low mass loss rates are necessary to re- 
produce the observed emission strength, i.e., the larger the 
opening angle, the lower the equatorial densities for a given 
emission strength. Also, the lower the density, the farther 
away from the star is the dust destruction surface and, there- 
fore, the lower the density in the dust condensation zone. The 
two effects combined have the consequence that models with 
A^op si 15° have densities below the critical density neces- 
sary to form dust (eq. lfT2l ). 

Models with small opening angles require high equatorial 
densities to reproduce the line emission an d, th erefore, can 
form dust. However, as we shall see in Sect. 13.21 in our mod- 
els dust is typically confined within a wedge about one open- 
ing angle away from the equator. This is because the den- 
sity drops quickly towards the pole and soon falls below the 
critical density for grain formation. If the opening angle is 
too small (A6'op < 5°), the dust will encompass too small a 
solid angle for the radiation that comes from the inner enve- 
lope (both stellar radiation and reprocessed radiation by the 

so deeply embedded in the C S out flow that differences in the emerging flux 
would be negligible (see Sect. l3.H . 
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Fig. 1. — Comparison between observations and the model SED for the 
three dust models used. The observations (red lines and symbols) are de- 
scribed in Paper III. Symbols: circles - UBVR photometry, squares - near-IR 
photometry from 2MASS, triangles - BASS data, large open circles - MSX 
data, large open squares - IRAS data, and solid lines - 2003 NIRIS data. Top: 
Dust model 1. Middle: Dust model 2. Bottom: Dust model 3. In all panels 
the dashed Hne represents the unattenuated stellar spectrum. Vertical hnes on 
the best-fit model SED are hydrogen emission lines (from left to right - H/3, 
Ha, and Br7). 

CS gas) and, as a consequence, will reprocess little radiation. 
This is a geometrical effect and is true even for very large dust 
opacities. Therefore, small opening angles can be ruled out on 
the grounds that they produce too small IR excess. 

Another trend that became evident is that dM{Q)/dVt cannot 
be larger than about 3 x 10"^ M0 yr"' sr"'. A large density in 
the poles implies that the polar regions would contribute sig- 
nificantly to the overall line emission. Since the polar material 
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Fig. 2. — Comparison between the model (black) and observed (red) line profiles. The observations and instruments used are described in Paper HI. Left: 
Best fit for the 2006, December 13 epoch, using dM(90°)/dn = 6 X 10"' Mq yr"' sr"' . Right: Best fit for fire 2006, December 27 epoch, using dM(90°)/dn = 
1 X 10"* Mq yr"' sr"'. 



has a large terminal speed, the polar emission would appear 
as a broad red emission wing, which is not observed. The po- 
lar emission begins to alter the emission profile significantly 
if mass loss rate exceeds the quoted value. 

This upper limit for dM{Qi)/dQ, has an important implica- 
tion on the possible range of inclination angles. Since the 
polar regions must be relatively optically thin, they cannot ac- 
count for the broad and deep absorption component observed 
in the Balmer lines (see, for instance. Fig. 3 of Paper III). 
Models with < / < 75° do not produce the necessary ab- 
sorption and, therefore, this range of inclination angles can be 
excluded. Besides this interval, our models can also exclude 
inclination angles too close to edge-on (/ w 90°). In all our 
models the radial optical depth along the equator is very large 
(r w 10^ in the Lyman continuum and k, 10^ in the Balmer 
continuum). As discussed in Paper III, some features of the 
stellar photosphere are observed, which, therefore, excludes 
edge-on viewing. 

Our best-fit SED, Ha, H/3, H7 and W5 fine profiles for the 
parameters listed in Table [T] are compared to the observations 
in Figs.[T]to[3] The CS gas density and velocity profiles along 
several latitudes are shown in Fig.|4] 

As said above, the opening angle and the equatorial mass 
loss rate are the primary parameters that control the strength 
of the line emission for a given /. Our best-fit value for the 
opening angle is 7°, which corresponds to m = 92. For the 
equatorial mass loss rate our best-fit value is in the range 
dM{9Q°)/dn = (1.4 — 1.6) x lO-** M© yr"' sr"'. The lower 
limit corresponds to the best-fit value for the line emission 
observed in 2006 December 13, and the upper limit repro- 



duces better the 2006 December 27 observations, which has 
substantially stronger emission (see Fig. |2]i. This result indi- 
cates that IRAS 00470-1-6429 presents significant changes in 
the mass loss rate in very short timescales (two weeks). 

Both m and dM{9Q°) / dO, are well constrained parameters 
in our models, but we note that reasonable fits of the emission 
strength could also be obtained with slightly larger opening 
angles (say, by one or two degrees) and slightly smaller (a 
few tens of per cent) dM{9Q°)/dVi. The opposite is also true, 
i.e., similar fits can be obtained by slightly smaller opening 
angles and slightly larger mass loss rates. 

The opening angle and the mass loss rate also control the 
SED shape. Their effects on the SED are complex, and de- 
pends significantly on the adopted /. Dust also plays a major 
role in defining the shape of the SED, not only in the IR but 
also on the visible spectral region provided that the line of 
sight of the observer crosses dusty material. For this reason, 
we postpone the discussion of th e model SED and its compar- 
ison to observations to Sect. 13.21 

Contrary to the equatorial mass loss rate, the polar mass 
loss rate is rather poorly constrained. As explained above, an 
upper limit of about dM(0)/dfl = 3x10"^ M© yr"' sr"' could 
be set on the basis of the absence of a broad emission com- 
ponent. Given the best-fit value of dM{9Q°)/dVl, above, this 
corresponds to a lower limit of Ai = 19, which means that 
the ratio between the equatorial and polar mass loss rates is 
at least as high as 20, but our models do not exclude much 
larger values of A 1 . The reason why our models do not con- 
strain well the polar densities is that the H I recombination 
lines are not suitable probes of low density regions. UV reso- 
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Fig . 3 . — Comparison between the observed H7 and US line profiles (2006, 
December 13th epoch, red) and models (black). 

nance lines of high ionization species are the ideal probes of 
those regions, but unfortunately no UV data are available for 
this object. 

It should be noted that even though the lower limit for the 
ratio between the equatorial and polar mass loss rates is 20, 
the ratio of the actual densities is much larger. As shown in 
Figure |4] this ratio is about 200 or larger, depending on the 
adopted value for A 1 . 

The shape of the absorption component of the H I lines is 
highly variable, with the blue absorption edge moving from 
« 400 to « 800 km s"' on very short timescales (see Fig. 3 
of Paper III). For this reason, it is difficult to reproduce the 
detailed shape of the absorption profile. The high variabil- 
ity of the observed line profiles suggests variable mass loss 
on short timescales. For instance, our models reproduce very 
well the H/3 absorption profile of 2006, December 27, but not 
the profile that was observed two weeks earlier. This may in- 
dicate that not only the mass loss vary but also the kinematic 
properties of the CS envelope change as well. 

The depth of the absorption component, contrary to its 
shape, is remarkably constant. It depends on two competing 
processes: the absorption of resonant radiation by H I atoms 
in a given excitation state and the emission filling-in, i.e., the 
generation of new photons both by line emission and bound- 
free and free-free continuous emission that partially fills in the 
absorption component. The filling-in by continuous emission 
is particularly important because of the very large densities 
close to the star (see FigurelDi. As shown in Fig. [5] our mod- 
els also fit reasonably well the depth of the absorption compo- 
nent for other lines of the Balmer series, which is not a trivial 



result. 

Another important parameter is the value of /3 that tells how 
fast the wind accelerates. We found that the best models for 
the line profiles are the ones with large values of /3 for all 
latitudes. Low values of (3 (say, /3 < 2) would result in a too 
broad emission component, due to the large velocities in the 
inner part of the wind, and also a in sharp absorption edge at 
high velocities. Both these features are not observed. Even 
though in our models it was possible to have different values 
of f3 for different latitudes, our best fit value was for /? = 4 and 
A3 = 0, i.e., a constant value of f3. Again, this is a somewhat 
uncertain result. For instance, a model with A3 = 1/3, for 
which /3{Q) = 3, is also a good model and cannot be ruled out 
from the present analysis. 

We end this section discussing the temperature and ioniza- 
tion structure of the CS gas. Fig.|5]shows the midplane tem- 
perature and the fractional number of H atoms in the ground 
level as a function of the distance from the star Clearly, 
the gas in hig hly non-isothermal, a resu lt already obtained by 
iDrewl (|I985[) and, more r ecentl y, by of ICarciofi & BiorkmanI 
(I2OO6I) and IZsargo et alj (l2008h . The temperature reaches a 
minimum of about 5 000 K at a distance of around 5 R and 
then steadily rises to an equilibrium temperature of about 
12000 K for r > 100 R. This rise in the temperature is prob- 
ably a result of our unrealistic chemical composition. The 
inclusion of the cooling terms by other elements will result i n 
a different temperature profile (see, for instance. lDrew|[T985l) . 

The gas is essentially ionized everywhere in the envelope. 
The maximum value for the relative fraction of H I and H II 
is never larger than 10%. Th is result agrees with the recent 
work by IZsargo et al.l (120081) . who studied the envelopes of 
sgB[e], and found that H is generally ionized, although this 
may depend on both the effective temperature of the central 
star and on the density scale of the wind. 

3.2. Circumstellar Dust 

Many processes are involved in defining the shape of the 
SED, and both gas and dust can play an important role de- 
pending on the wavelength range considered. In the visi- 
ble region, the dominant processes are continuum absorption 
by both gas and dust, and emission by continuum processes 
in the gas (bound-free and free-free emission). The SED is 
controlled by both gas and dust emission at wavelengths of 
< 2 /im in the near-IR, but the dust becomes the major player 
at longer wavelengths 

For this initial study of the dusty content of 
IRAS 00470+6429 three different models for the CS dust 
were investigated. For all models the grain size distribution 
has the form of eq. ( fTST l with n = -3.5, and each model ex- 
plores different values for the minimum and maximum sizes 
of the distribution. The first model studied correspond s to the 
standard MRN model for interstellar dust grains (Mat his et alJ 
Il977h . for which a^^ = 0.05 and a^ax = 0.25 /im. The second 
model employs a much larger upper limit for the distribution 
(flmax =10 /im), whereas the third model has very large grains 
only (amin = 1 and Umax = 50 /im). The three dust models are 
listed in Table [T] 

At this point, it is useful to consider what controls the shape 
and level of the IR excess arising from dust thermal emission. 
The shape of the emitted spectrum is controlled mainly by 
the dust absorption coefficient in the IR domain, which is a 
function of chemical composition and size of the grain. Other 
factors controlling the shape are the temperature and density 
structure. The level of the IR excess depends on the fraction 
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Fig. 4. — Top: gas density vs. distance to tlie star, for several colatitudes, 
as indicated. 6 = 90° corresponds to the equator. Bottom: radial velocity vs. 
distance. 

of the radiation coming from the central regions that is re- 
processed by the dust grains. Since in our model the dust is 
confined to the dense equatorial outflow, the amount of repro- 
cessing will be a function of the disk opening angle, A^op, 
and the dust optical depth, which is set by the gas-do-dust 
ratio and pdust (the bulk density of the grain material). 

The role of t] in defining the dust condensation site deserves 
a more careful examination. Recall that the dust condensation 
site is defined by two conditions: the density of the grow- 
ing material (Si) must satisfy eq. (fT2] i. in which 77 is a scaling 
factor, and the dust equilibrium temperature must be smaller 
than a given value that corresponds to the temperature above 
which the dust grains are destroyed (Sect. l2!2] i. In Fig.|6] we 
plot the equilibrium temperature of the smallest grains con- 
sidered by us (fl = 0.05 jim). As we see below, the smallest 
grains are always hotter than the larger grains, and, therefore, 
are used as the tracers of the dust formation site. The plot 
represents a cut through the envelope such that the direction 
of the symmetry axis is vertical and the direction of the mid- 
plane is horizontal. The latitude-dependent inner dust radius 
represents our solution for the dust destruction surface. Near 
the midplane, where the densities are larger and the material is 
more shielded from the central radiation, the dust destruction 
surface is closer to the star, and this surface moves outward 
for higher latitudes. 

The solid white lines in Fig. |6]mark the points of the wind 
for which the right side of eq. ( fT2] i is equal to the left side, 
and it encloses the region around the equator which have den- 
sities larger than the critical density for grain formation (to 
draw this line we used 7] = 0.02, see below). According to our 




Fig. 5. — Top: electron temperature along the midplane vs. distance to 
the star. Bottom: ratio between the number density of H atoms in the ground 
level, »i, to the total number of H atoms (im) and protons (rip). Shown are 
the results for dM(90°)/dn = 6 X IQ-' Mq yr"' sr"'. 
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Fig. 7. — Dust equilibrium temperature as a function of distance to the 
star. Top: dust model 1. Bottom: dust model 2. In each panel we show the 
temperature for the smallest (solid line) and largest (dashed line) grains of the 
distribution. 

two assumptions, the site for dust formation is the region of 
the wind both beyond the destruction surface and inside the 
region enclosed by the solid white lines. The dust condensa- 
tion site has an angular size of 2A6'(just, the value of which is 
controlled by the parameter r/: the larger 7] the smaller A6'dust- 

It is important to note that the site for dust formation does 
not necessarily coincides with the site where dust actually is, 
because the dust formed in the innermost region of the con- 
densation site will be carried outwards along the wind stream- 
lines. Therefore, dust that was formed at high latitudes will 
be dragged to regions where the density is much too low to 
form dust. This mechanism allows dust to exist at high lati- 
tudes in the outer envelope, where the density is lower than 
the critical value. Therefore, in our models dust is confined 
within a wedge with a (half) opening angle of A6'dust- A^dust 
is a function of both A^op and 7]. 

Figure [Ucompares the best-fit SED obtained for each of the 
three dust models tested. We did not attempt to fit the IRAS 
fluxes, but focus instead on the BASS and MSX data, since 
the large aperture of IRAS has probably lead to some degree 
of contamination by the interstellar cirrus. Also, we recall, 
from the data description in Paper III, that the NIRIS data 
have large errors of up to 20% in the photometric calibration. 

The MRN model (model 1, top panel of Fig.[T]i does not fit 
well the data. The integrated IR flux is in rough agreement 
with the observations, which means that the model does re- 
process the correct fraction of the bolometric stellar luminos- 
ity. It reproduces the SED reasonably well in the 2 — 3 fim 
range, but fails to reproduce the slope of the SED for larger 



wavelengths. The model with a large upper limit for the grain 
size (model 2, middle panel of Fig. ^ fits both the level and 
shape of the SED in the 2 — 7 fim range. The reason why 
model 2 fits better the shape of the near-IR SED than model 
1 resides in an interesting physical effect: differentl y sized 
grain s have different equilibrium temperatures (Carciofiet al.l 
120041) . Figure |7] shows the temperature of the smallest and 
largest grains of dust models 1 and 2. For Model 1, the span 
in grain sizes is not very large and the difference between the 
coolest and hottest grains is, at most, 300 K. For Model 2, 
however, we find that the very large grains are much cooler 
then the small grains. This explains the difference in the slope 
of the SED: model 2, having on average a lower dust equi- 
librium temperature, have proportionally larger emission at 
longer wavelengths. 

The 10 fim region of the spectrum presents a big challenge 
for models 1 and 2, since for both the 9.7 /^jm silicate emis- 
sion feature is quite strong, in sharp contrast to the observa- 
tions. This situation is much improved for the model with 
only large grains (model 3, bottom panel of Fig. [T]). This 
model reproduces reasonably well the shape of the thermal 
IR emission in the entire 2 — 13 /im range. As discussed 
above, the shape of the IR thermal continuum is controlled 
mainly by the dust absorption opacity, which is a function of 
grain size and composition. Under our assumption that the 
grains in IRAS 00470+6429 are oxygen-based, the observed 
IR excess seems to indicate that large grains {a> I /im) are 
the dominant constituents of the CS dust. 

The large difference between the equilibrium temperatures 
of small and large grains have another consequence for the 
models. Recall that one of the conditions for determining the 
dust condensation site is that the temperature of the hottest 
grains cannot exceed rdestruction- For models 1 and 2 this con- 
dition is satisfied only for r > 540 R, but for model 3 the dust 
destruction surface is much closer to the star, r>19R. 

Another important difference between models with small 
and large grains lies in the bulk density of the grain mate- 
rial (pdust)- At a given location of the envelope, a constant 
fraction of the gas is converted into dust. Since the mass of 
the grain increases as a^ and its cross section as ^ a^, we 
had to decrease the bulk density of the material so that the 
dust optical depth remains approximately constant. There- 
fore, our modeling for the IR excess of IRAS 00470+6429 
suggests that the dust grains are both large {a > I /im) and 
fluffy (pdust -- 0. 1 g cm"^). 

Finally, the three models also reproduce well the visible 
SED (Fig.[T]i, but require different amounts of interstellar red- 
dening: for the models with small grains, £'(B-y)is = 1.1, and 
for the model with large grains, E{B -V)is = 1.2. Recall that 
E(B - y)is is a free parameter in the modeling and was set by 
fitting the visible SED. The different values of ^(B-y)^ can 
be easily understood by considering that the CS extinction of 
model 3 is essentially gray in the visible and for this reason its 
intrinsic CS reddening is smaller than the models with small 
grains. 

The extent to which the CS disk modifies the stellar radi- 
ation can be assessed by comparing the model SED with the 
unattenuated stellar radiation (dashed line in Fig. [TJ. Since 
the star is viewed close to edge on, this causes the extinction 
of UV and visible radiation, that is reprocessed into longer 
wavelengths. Most of the reprocessed radiation escapes in the 
polar direction, due to the fact that the vertical optical depth 
is much smaller than the radial optical depth, but part of it 
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escapes radially. The [/-band flux is particularly important 
for the SED fitting. Our models predict that when the object 
is viewed pole-on, the large bound-free excess produced by 
the inner envelope completely fills in the photospheric Balmer 
jump, thus producing a power-law SED up to the UV region. 
The fact that the Balmer jump is present is another evidence 
in support of larger inclination angles (Sect. 13.1b . 

4. DISCUSSION 

At this point it is useful to recall the assumptions of the 
model presented here. For this exploratory study it was as- 
sumed a bimodal wind composed of a slowly outflowing disk 
and a fast polar wind, and an ad hoc density and velocity law 
was used to explore different characteristics of this bimodal 
wind model, such as the ratio between the mass loss rate at 
the pole and the equator, and the opening angle of the dense 
equatorial outflow. Here we intentionally limited the number 
of free parameters and focussed our analysis on the CS region, 
aiming at studying its gaseous and dusty phases. The use 
of the radiative transfer code HDUST allowed us to perform 
a fully consistent solution of the radiative transfer problem. 
Since the results are not marred by unknown radiative transfer 
effects, the comparisons between the models and observations 
represent genuine tests of the physical assumptions. 

In our models the line emission comes mostly from the 
dense outflowing disk, which has an opening angle of A6'op = 
7°, whereas the strong line absorption component originates 
in the upper layers of the disk. The model line profiles re- 
produce the emission component of the H I lines, as well as 
important features of the absorption component, such as its 
depth, reasonably well. 

The total stellar mass loss can be calculated by integrating 
eq. (|2]i over all solid angles 



M = 2ti 



^(0)sin(0M0, 
ail 



(14) 



Using the best-fit values from Table [T] we obtain M = (2.5 
— 2.9) X 10"^ Mq yr"'. The lower limit corresponds to the 
fit for the 2006, December 13 epoch, whereas the upper limit 
describes the 2006, December 27 epoch. The total mass loss 
rates derived here are one to two o rders of magnit ude lower 
than what is estimated for sgB[e] dZickgraf) 120061) . On the 
other hand, it is at least two orders of magnitude larger than 
the mass-loss rate of main-seque nce B stars, which is of the 
order of 10"'° — 10"'' Mq yr"' (iGafliier et al.lfMl . This 
result, if corroborated by future studies, leads to the question 
of how mass loss is driven in FS CMa objects. 

An inclination angle / k, 84° was obtained from the model- 
ing. Since the system is viewed nearly edge-on and the emis- 
sion component of the line profile is narrow, the presence of 
a large rotational component in the velocity field can be ex- 
cluded, thereby justifying the pure radial expansion assumed 
here. It follows that a rotationally supported Keplerian disk, 
which i s one of the scenarios proposed for the disks around 
sgB[e] (lPortedl2003h . can be excluded for IRAS 00470H-6429. 

According to our assumptions for grain formation, the dust 
around IRAS 00470H-6429 is contained within a wedge of 
opening angle A^dust ~ 10°. This opening angle was found by 
fitting the IR excess, which is a direct measure of the amount 
of radiation that is reprocessed by the CS dust. This value 
of A^dust implies that the scaling factor -q in eq. ( fT2] i is about 
0.02, therefore much less than unity. The physical implication 
of this value of rj is that the critical density for grain formation 



derived by iGail & Sedknavil ( 119881) might be overestimated, 
i.e., grains might be able to form on shorter timescales. 

The models clearly indicate a prevalence of large grains in 
the dust around IRAS 00470H-6429, for the reason that model 
3 was the only model to satisfactorily fit the IR excess and 
to explain the absence of a conspicuous 9.7 /im silicate emis- 
sion feature, which is to be expected for oxygen-based dust 
grains. This result seems to be in contradiction with the usual 
view on grain formation, according to which large grains form 
from the coalescence of smaller grains. This apparent con- 
tradiction can be explained by assuming that the medium 
around IRAS 00470H-6429 is clumpy. In such a medium, 
dense clumps are the ideal sites for the grain formation and 
growth, because they provide the necessary shielding from 
the surrounding radiation. As the clumps expand, follow- 
ing the wind streamlines, they become optically thiner, ex- 
posing the grains to more radiation. Because small grains 
are hotter than large grains, they will have a larger proba- 
bility to be destroyed . Clumps are thought to exist in the 
wind of hot stars (e.g. lKudritzki & Pulsll200(3l) . and this sim- 
ple mechanism of selective grain sublimation may be able to 
explain the observed trend of a prevalence of large grains in 
IRAS 00470H-6429. A careful investigation of this hypothesis 
will be left for a future study. 

Another implication from our modeling is that there is a 
significant CS extinction in the visible and UV regions. This 
follows from the almost edge-on orientation of the equato- 
rial plane of the CS envelope. As a result, the object looks 
fainter than without the envelope, which reduces its estimated 
distance. The distance that follows from our modeling is 
d ^\.\ kpc, depending on the dust model. However, the large 
E{B-V)i^ that came from the modeling of the SED in the vis- 
ible range, 1.1 — 1.2 mag, is more appropriate for a larger 
distance (~ 2 kpc). A similar distance estimate was also sug- 
gested in Paper III (2.0±0.3 kpc), but in that work the CS 
reddening was not accounted for. This apparent discrepancy 
between distance and E{B-V)is may be reconcile d by in- 
creasing the stellar luminosity. As mentioned in Sect. 12.51 the 
stellar parameters are only loosely constrained by our model- 
ing, and a higher luminosity (and therefore a larger distance 
to the star) is certainly plausible. 

A critical parameter in our modeling is the inclination an- 
gle. Better constraints on this parameter will be put after po- 
larization measurements become available. Modeling the po- 
larization along with the SED and spectral line profiles can 
be done with HDUST and this will allow us to constrain much 
better the geometrical properties of the CS envelope, includ- 
ing the inclination angle. 

5. CONCLUSIONS 

In this paper we adopted an ad hoc model for 
IRAS 00470H-6429 which is physically motivated by the cur- 
rent knowledge of the CS envelopes of sgB[e]. The model 
consists of a 6 Rq central star with Teff = 20000 K and L = 
5 100 Lq, surrounded by a bimodal CS envelope composed 
of a dense, slowly outflowing disk-like wind and a fast polar 
wind. The radiative transfer calculations were performed with 
the computer code HDUST . 

In the model the underlying physical reason for the bimodal 
envelope is that the mass loss is, somehow, enhanced around 
the equator. From the analysis of the observed SED and H I 
line profiles, it was found that the equatorial outflow is at least 
200 times denser than the fast polar wind and about 10 times 
slower. In addition, the opening angle of this slowly outflow- 
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ing disk (defined as the co-latitude for which the mass loss 
rate drops to half of its value at the equator) is about 7°. Our 
analysis firmly exclude both very small (< 5°) and very large 
(> 15°) opening angles. 

We determined that the integrated stellar mass loss rate is 
M = 2.5 — 2.9 X 10"^ Mq yr~'. The two extremes correspond 
to fits to observations made in December 2006, but separated 
by an interval of two weeks. This gives a quantitative measure 
of how the mass loss varies in short timescales. The above 
value of M, while much smaller than that of sgB[e], which 
is of the order of 10"'' — lO"*" M© yr"', is at least 100 times 
larger than that of main-sequence stars of spectral type B. 

We adopted a prescription for dust formation based on two 
complementary criteria. In order to form dust at a given point 
in the CS envelope, the equilibrium temperature of the dust 
grains must be smaller than a given grain destruction tem- 
perature, Tdestniction, and the gas density must be larger than a 
critical value. Several important radiative transfer effects are 
considered in our calculations, such as the shielding of the 
dust by the optically thick inner CS material and the fact that 
differently sized grains have different equilibrium tempera- 
tures. 

To investigate the properties of the dusty content of 
IRAS 00470+6429 we studied dust models with three differ- 
ent grain size distributions, one with the standard MRN dis- 
tribution (fl = 0.05 — 0.25 /im), one with both small and large 
grains (a = 0.05 — 10 /im), and one with only large grains 
(a = 1 — 50 /im). The dust was assumed to be oxygen-based 
(silicates), an assumption supported by evolutionary argu- 



ments and our Spitzer observations of FS CMa stars. 

The only model capable of reproducing the observed IR 
excess in the entire 2 — 13 /im range was the model with 
only large grains (omin = 1 /im), because the presence of small 
grains always results in a strong 9.7 /im silicate emission, 
which is not observed. One consequence of the prevalence of 
large grains around IRAS 00470+6429 is that the bulk density 
of the grain material must be very small (pdust '^ 0. 1 g cm""*). 
Therefore, the observed shape of the IR excess seems to firmly 
indicate that the CS dust grains of IRAS 00470+6429 are both 
very large and fluffy. 
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